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Detection of a new, low-brightness supernova remnant possibly 
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Abstract. We report on the discovery of a shell- type supernova remnant in the southern sky. It is a large (~ 8° x 8° ) , 
low-brightness source with a nonthermal radio spectrum, which requires background filtering to isolate it from 
the diff'use background emission of the Galaxy. Three 3EG 7-ray sources are spatially correlated with the radio 
structure. We have made 21-cm line observations of the region and found that two of these sources are coincident 
with HI clouds. We propose that the 7-ray emission is the result of hadronic interactions between high-energy 
protons locally accelerated at the remnant shock front and atomic nuclei in the ambient clouds. 
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1. Introduction 

The origin of cosmic rays (CRs) is a long-standing problem 
in contemporary astrophysics. It has been thought since 
the 1960s, on the basis of energetic considerations, that 
protons are accelerated up to energies of ~ 10^^ eV in su- 
pernova remnants (SNRs) (Ginzburg & Syrovatskii 1964). 
When first-order Fermi acceleration at shock fronts was 
established as an efficient mechanism for CR production 
(e.g. Bell 1978), SNRs became the most natural candi- 
dates for the accelerating agents of the bulk of galactic 
CRs. Multifrequency radio observations clearly show that 
SNRs accommodate a population of relativistic electrons 
with the expected energy distribution. However, there is 
at present not such a conclusive observational evidence for 
shock acceleration of protons and ions near SNRs. 

High-energy 7-ray astronomy is the most appropriate 
tool to probe SNRs in the search for an observable signa- 
ture of proton acceleration. The interactions between rel- 
ativistic protons and ambient nuclei can produce neutral 
pions (tt^) which decay to yield 7-ray emission that might 
be detected by instruments like the Energetic Gamma Ray 
Experiment Telescope (EGRET) on board the Compton 
satellite. Possible associations of 7-ray sources with SNRs 
have been claimed by several authors in the past (Pollock 
1985, Combi & Romero 1995, Esposito et al. 1996, and 
Combi et al. 1998a). However, a leptonic origin of the 7- 
ray emission (e.g. Pohl 1996) cannot be ruled out in most 
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of these cases, especially when relatively young remnants 
are considered. 

In this paper we present results of a new observa- 
tional study designed to clarify the problem of the origin 
of galactic CRs. Aharonian et al. (1994) and Aharonian 
& Atoyan (1996) have pointed out that the presence of 
nearby clouds can increase the probability of detecting 
the nucleonic CR component locally accelerated in a SNR 
because the clouds act as targets for the relativistic pro- 
tons, hence enhancing tt" production and the resultant 7- 
ray luminosity. A large remnant in a cloudy medium may 
have several discrete 7-ray sources associated with it. The 
study of the spectra and fluxes of these 7-ray sources could 
provide information about the cosmic ray production at 
several different sites where the clouds are overtaken by 
the accelerated protons. 

Some of the 170 unidentified 3EG 7-ray sources 
(Hartman et al. 1999) appear to occur in small clusters 
subtending less than 6° in the sky at high galactic lat- 
itudes. We have searched the 408 MHz all-sky survey 
(Haslam et al. 1981) for radio continuum emission to- 
wards the directions of those four clusters of high- latitude 
sources that are visible from the southern hemisphere. 
The EGRET confidence contours for the cluster consist- 
ing of 3EG J1834-2803, 3EG J1847-3219, and 3EG J1850- 
2652 correlate spatially with an extended shell-type ra- 
dio source, which is masked by steep gradients in the 
diffuse galactic radiation. We have subsequently per- 
formed Hl-line observations towards this large area in the 
Capricornus region, and have detected hydrogen clouds at 
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Fig. 1. Background filtered radio emission at 408 MHz of 
the region surrounding the three 7-ray sources. Contours 
are labelled in steps of 1 K in brightness temperature, 
starting at 1.2 K. The superposed grey-scaled levels rep- 
resent the 99%, 95%, 68%, and 50% statistical probability 
that a 7-ray source lies within each contour according to 
the EGRET catalog (Hartman et al. 1999). HPBW ~ 5l' 

the positions of two of the 3EG sources. Additional radio 
continuum data was provided by the Rhodes/HartRAO 
2326 MHz survey (Jonas et al. 1998). Having radio im- 
ages at two widely-spaced frequencies allows the spatial 
distribution of the spectral index to be computed, hence 
establishing the nonthermal nature of the emission and 
revealing the energy distribution of the electronic compo- 
nent of the CRs at the SNR. In what follows we describe 
our data analysis procedures, new observations, and dis- 
cuss the results of our study. 

2. Data analysis and new observations 

The detection of weak, extended radio sources is a diffi- 
cult task due to the contamination produced by the diffuse 
background radiation in the radio images. Faint sources 
with diameters of several degrees can only be identified 
by using effective spatial filtering techniques applied to 
good-quality, single-dish data. In this study we have used 
data at 408- MHz and 2326-MHz from the radio surveys by 
Haslam et al. (1981) and Jonas et al. (1998), respectively. 
The angular resolutions (HPBW) of these surveys are 51' 
(408 MHz) and 20' (2326 MHz). We have extracted maps 
of 12° X 12° around the mean position of the 7-ray sources 
and have then used the background filtering method orig- 
inally developed by Sofue & Reich (1979) in order to em- 
phasize the fine radio structure superposed on the diffuse 
galactic radio emission. This iterative procedure results in 
the elimination of diffuse sources with angular scales larger 
than a specified filtering beam (see Combi et al. 1998b for 
more details). This technique is very versatile, can be used 
at any galactic latitude, and is well-proven for a range of 
angular scales. It was recently used to discover the radio 
counterpart of the nearby SNR RX J0852.0-4622, which 
is contaminated by the bright, extended emission of the 
Vela region (Combi et al. 1999). 

When the 12° x 12°-maps are filtered, a steep gradient 
of emission shows up towards the galactic plane. In order 
to remove this residual background, we filtered a larger 
region of 20° x 24°, which includes the galactic plane. Then 
the smaller maps were extracted again from the filtered, 
larger map. 

Figures 1 and 2 show the background filtered radio 
emission of the region under study at 408- and 2326- 
MHz, respectively. A filtering beam of 240' x 240' was 
applied to both images (i.e. all emission structures with 
angular scales larger than 4° have been suppressed in 
these images). The probability confidence contours of the 
3EG sources have been superposed in the maps. A large 



Fig. 2. Idem Fig. 1, but at 2.326 GHz. Brightness tem- 
perature contours are in steps of 0.03 K, starting at 0.03 
K. HPBW ~ 20' 

(^8° x 8°), shell- type source can be clearly seen at both 
radio frequencies. In order to determine the spectral index 
distribution of the shell emission we have convolved and 
re-tabulated the 2326-MHz data to the same beam and 
tabular interval as the 408-MHz map. Following Combi 
et al. (1998b), the spectral indices and the correspond- 
ing errors were estimated from these maps. Uncertainties 
from the filtering method also should contribute to the 
errors, but are difficult to estimate. We think that a con- 
servative assumption here is to assume that the real errors 
double the formal ones. The spectral index distribution, 
which is shown in Figure 3, is consistent with indepen- 
dent estimates obtained using a different filtering method 
developed by Jonas (1999). 

Aharonian & Atoyan (1996) emphasize that the ex- 
istence of a potential CR accelerator, such as a SNR, is 
not a sufficient condition for the generation of detectable 
levels of 7-ray emission. Target nuclei are also required, 
so we have performed HI observations towards the area of 
interest in order to search for ambient atomic gas in the lo- 
cal neighbourhood of the SNR. The lAR 30-m single dish 
telescope at Villa Elisa, Argentina was used for these ob- 
servations. The receiver was a helium-cooled HEMT am- 
plifier with a 1008-channel autocorrelator at the backend, 
which has been used for similar purposes by Combi et al. 
(1998a, b). System parameters and additional details of 
the observational technique can be found in these papers. 

The HI line was observed in hybrid total power mode 
and the sky was sampled on a 0.5° rectangular grid. Each 
grid position was observed for 36 s with a velocity resolu- 
tion of ~ 1 km s~^ and a coverage of ±450 km s~^. A set 
of HI brightness temperature maps (ATrms ~ 0.2 K) were 
built for the velocity interval (—30 km s~^, -|-50 km s~^), 
which is appropriate for this study. We adopted the fol- 
lowing criteria to identify standard clouds in this velocity 
range: the maps should show a closed, simply connected 
structure (i.e. no shells or filaments were considered), with 
a clear maximum in brightness temperature, and an ex- 
tension in velocity of at least 6 km s~^. There are six 
structures that satisfy this criterion in the region of inter- 
est. The two clouds, CI centred at {I, b) ~ (10°, -12°) 
and C2 centred at (/, b) ^ (5.5°, —8°), are the stronger 
ones, with brightness temperatures of 45 K, which are 
detected at a confidence level of 225cr. Figure 4 shows 
a series of HI brightness temperature maps where these 
clouds can be clearly appreciated. In Figure 5 we present 
the integrated column density map for the relevant veloc- 
ity interval with the 3EG sources superposed. Two of the 
sources (3EG J1834-2803 and 3EG J1850-2652) are posi- 
tionally coincident with the clouds CI and C2. Standard 
galactic rotation models (Fich et al. 1989) indicate that 
the remaining clouds, not shown in Figure 4, are located 
at kinematic distances that are far too distant to be phys- 
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Fig. 3. Spectral index distribution computed between 
408 MHz and 2.326 GHz for the region shown in Figs. 1 
and 2. The EGRET probabiUty contours are superposed. 

Table 1. Measured properties of the new SNR in 

Capricormis 



Property 


Value 


Galactic coordinates (center) 


(+6.5°, -12.0°) 


Angular size (deg xdeg) 


8.0 X 8.0 


Flux density (408 MHz) 


~ 180 Jy 


Flux density (2326 MHz) 


~ 55 Jy 


Average spectral index 


-0.68 ±0.15^ 


^Formal error. 



ically related to the continuum radio source, which due 
to its largo size is constrained to be nearby. It is also in- 
teresting to notice that there is a local hole in the HI 
distribution at the approximate position of the center of 
the SNR. 

3. Main results 

The newly-identified extended shell-like radio source 
in the Capricornus region is centred at {l,b) « 
(+6.5°, —12.0°). It has an integrated flux density of 
^ 180 Jy at 408 MHz, once the contribution from back- 
ground point sources has been removed. This extragalac- 
tic contribution is dominated by the sources PKS 1839-34, 
PKS 1848-338, and PKS 1853-300. Flux errors arc diffi- 
cult to estimate due to uncertainties related to the filtering 
technique, in such a way that the quoted flux should be 
considered only as a rough value. 

The identification of the extended source as a new SNR 
is not only supported by its large size and shell-like mor- 
phology, but also by the spectral indices of the radio emis- 
sion (see Fig. 3) which are clearly nonthermal with an 
average value of a « —0.68 ± 0.15, somehow steeper than 
those observed in typical of shell-type remnants (Green 
2000). Notice, however, that the quoted errors are purely 
formal and, as it was previously mentioned, the real ones 
could be significantly larger (perhaps around 0.3). The 
measured properties of this new Capricornus SNR are 
summarized in Table 1. 

The distance to the remnant is uncertain, although 
its large size and high latitude suggest that it is nearby. A 
rough distance estimate of ~ 470 pc can be obtained us- 
ing the S(Z))-dependence for low-brightness shell-like rem- 
nants proposed by Allakhverdigev et al. (1986), although 
this method alone is too dubious as to provide a reliable 
result. An independent upper limit of 700 pc is provided 
by the relatively high value of |6|, if the source is assumed 
to be located within a galactic disk with a half-liciglit of 
150 pc. In addition, galactic rotation curve models (e.g. 
Fich et al. 1989) impose an upper bound of ~ 500 pc on 



Fig. 4. HI brightness temperature channel maps (contour 
labels in K) obtained for the velocity range —3.0 km s~^ 
to +7.0 km s~^. The supernova remnant boundary and 
the EGRET probability contours are superposed. 

the distance, assuming that the HI clouds Gl and C2 are 
physically related to the remnant. 

A radius of ^ 30 pc is obtained for the SNR if a dis- 
tance, not unreasonable, of 470 pc is adopted, indicating 
that it is probably at the end of its adiabatic phase. Using 
standard Sedov solutions for an initial energy release of 
^ 10^^ erg, we derive an age of ~ 16, 500 yr, assuming 
an average density of n ~ 0.05 cm~^ for the intercloud 
medium at high galactic latitudes. 

In Figure 4 we show the HI clouds Gl and G2 using 
an ensemble of channel maps for the velocity interval —3 
to +7 km s^^, but one of the clouds is visible even be- 
yond these velocity limits. Gl can be traced from —8 to 
+2 km s~^ whereas G2 extends from —1 to +5 km s~^. 
The cloud column densities were evaluated by integrating 
over these velocity intervals. 

The integrated HI column density maps were used to 
obtain estimates for the total masses of the individual 
clouds. Using the adopted distance of 470 pc we esti- 
mate masses of ^ 1400 Mq and ~ 1200 Mq for Gl and 
G2, respectively. The aggregate systematic error in these 
computed masses may be as large at 50 %. Dominant 
sources of error are uncertainties in the adopted distance, 
imperfect background subtraction, and possible existence 
of other (yet undetected) molecular species. 

4. Discussion 

The three 7-ray sources superposed with the SNR are clas- 
sified as being possibly extended in the 3EG catalog, and 
none of them show variability. Variability is not expected 
if the 7-ray emission is related to the radio remnant. 

The sources 3EG J1850-2652 and 3EG J1834-2803 are 
spatially coincident with the HI clouds Gl and G2, respec- 
tively. We suggest that the collective spatial coincidence 
of the 7-ray sources, HI clouds and the radio remnant pro- 
vides evidence for the production of 7-rays by 7r°-decay, 
where the pions result from proton-proton collisions at the 
SNR/ambient cloud interfaces. 

In order to make quantitative estimates of the proba- 
bility of chance association of the 7-ray sources in clusters, 
we have adapted the numerical code developed by Romero 
et al. (1999) and used by Romero, Benaglia, & Torres 
(1999) to study the positional association of unidentified 
EGRET sources with various populations of galactic ob- 
jects. The code calculates angular distances between dif- 
ferent kinds of celestial objects contained in selected cat- 
alogues, and establishes the level of positional correlation 
between them. Monte-Garlo simulations using large num- 
bers of synthetic populations are then performed in order 
to determine the probabilities of pure chance associations. 
When generating synthetic populations of 7-ray sources 
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the distribution in galactic latitude is constrained to be 
the same as that for the actual 3EG sources. This is nec- 
essary in order to obtain reliable results since the distri- 
bution of the 3EG sources is non-isotropic, with a strong 
concentration towards the galactic plane. The reader is re- 
ferred to the paper by Romero, Benaglia, & Torres (1999) 
for further details of the simulation code. 

Using this code we determined that there are 54 
unidentified 3EG sources (out of 116) that cluster within 
6° of each other at high latitudes (|6| > 5°). When 1500 
random populations of sources are generated (a larger 
number does not significantly change the results), the 
average number of sources in 6° clusters is found to 
be 37.6 ±4.7. This translates to a probability as low as 
1.9 X 10"'' that all high-latitude 6° clusters are chance 
events. For |&| < 5° there are 47 sources out of 54 in 6° 
clusters, whereas the simulations predict 34.6 ± 3.4 chance 
associations. The probability that all of the observed low- 
latitude clusters are chance events is 1.5 x 10^^, which is 
similar to the value obtained for |6| > 5°. These simula- 
tions show that several apparent clusters of unidentified 
7-ray sources are probably groups of physically related 
sources. 

We shall now estimate the probability of finding a SNR 
by chance in the direction of a southern cluster of 7-ray 
sources. In order to do this, we need to know the chance 
probability of finding a shell-type, relatively large SNR 
overlapping with a circle of the cluster size in an arbi- 
trary direction at intermediate (say 5° < |6| < 20°) galac- 
tic latitudes. There is just one of such SNRs in Green's 
(2000) catalog. Of course, many large and yet undetected 
remnants of low surface brightness should exist under the 
diffuse non-thermal galactic emission. A large-scale study 
with a systematic subtraction of the galactic background 
component of the southern radio sky at 2.326 GHz by 
Jonas (1999) has recently shown that the number of these 
uncataloged SNRs in the relevant latitudes is about 17. 
Since we have used data from the same Rhodes/HartRAO 
SKYMAP survey, we shall adopt this figure as a conve- 
nient estimate of the total number of remnants expected 
in the region of interest here. Then, the Poisson probabil- 
ity of finding one SNR overlapping with a cluster of 3 or 
more unidentified 3EG sources is ~ 1.4 x 10"^. 

The probability of chance association is further con- 
strained if we consider that there are two clouds appar- 
ently coincident with two of the three 7-ray sources in the 
cluster for which the new SNR was found. The probability 
p{n) that the line of sight along a distance r will intersect 
with n clouds is (Spitzer 1998): 



, , (kr)^e 

P{n) = j— 

ni 



where k is the average number of clouds per kpq^. In the 
direction of the galactic plane fc ~ 6.2 for standard clouds 
(Spitzer 1998). We found no published estimates of k for 
relatively high latitudes, so we adopt the galactic plane 
value in order to obtain an upper limit to the probabil- 
ity. Using this equation and a distance range equal to the 
estimated diameter of the SNR (60 pc), we obtain that 
the chance probability of finding one or more clouds in 
each of two different lines of sight is 0.097 (about 10 %). 
Consequently, the probability that the observed configura- 
tion of SNR, 7-ray sources and HI clouds be the sole effect 
of chance superposition is about 1.3 x 10^"^. This proba- 
bility, although not overwhelming, suggests that the hy- 
pothesis of a multiple 7-ray source that could be causally 
related to the SNR deserves serious study in the light of 
its important potential implications. 

There are mechanisms other than p — p induced tt*' 
decay that may contribute to the 7-ray emission from 
SNRs. Inverse Compton (IC) scattering of relativistic elec- 
trons can boost ambient low-energy photon fields (such 
as the cosmic microwave background radiation, the dif- 
fuse galactic infrared/optical radiation, and the remnant's 
own photon fields) to 7-ray energies (e.g. Pohl 1996, 
Mastichiadis 1996). Brcmsstrahlung emission from elec- 
trons in the shocked gas can also yield significant 7-ray 
luminosities (e.g. Blandford & Cowie 1982). The relative 
contribution of each of these processes to the overall 7-ray 
flux from the remnant depends on factors such as the age 
of the remnant, the injected particle spectra, the ambient 
gas density, and the ratio of electrons to protons at high 
energies (e.g. Sturner et al. 1997, Gaisser et al. 1998). 

In large and old remnants, such as the Capricornus 
remnant presented here, the IC boosting of diffuse back- 
ground photons should dominate that of local radiation 
fields (Sturner et al. 1997). At high energies electrons cool 
more efficiently than protons by synchrotron and IC losses, 
therefore the electron-to-proton ratio in the source is ex- 
pected to be ^ 1, i.e. it would be similar to what is ob- 
served in the CRs near the Earth (see Aharonian et al. 
1994, and especially Biermann & Strom 1993 for addi- 
tional arguments). IC radiation is therefore not expected 
to be dominant. Two factors indicate that the EGRET 
sources considered in this study are not predominantly 
due to IC scattering: (a) the detectable 7-ray emission is 
located towards regions where the observed radio emission 
(and hence the energy flux of relativistic electrons) is not 
particularly strong, and (b) two of the 7-ray sources occur 
towards HI clouds. If the HI clouds are sufficiently dense 
and the relativistic electron flux is high enough, then the 
observed 7-ray flux could be produced by bremsstrahlung 
interactions between the electrons and the ionized atoms 
in the compressed post-shock regions. This mechanism 
can, however, be discounted for our sources because both 

^ This is actually an oversimplification because it ignores 
the fact that clouds form a power-law distribution of prop- 
erties. The results, however, are not significantly changed in 
the present case, where massive clouds are considered. 
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Fig. 5. A map of integrated column density of HI (veloc- 
ity interval from —3.0 km to +5.0 km s^^) for the re- 
gion of interest. Label units are 10^^ atoms cm^^. EGRET 
probability contours are superposed. 

the cloud densities {n ~ 10 cm"'^, inferred from the HI 
data) and the electron flux (estimated from the radio 
synchrotron emission) are far too low. The electron flux 
necessary to produce the observed 7-ray luminosity by 
bremsstrahlung mechanism would result in a synchrotron 
radio flux of about 100 Jy at 1 GHz for each cloud. This ra- 
dio flux is two orders of magnitude greater than the value 
determined from the radio data for the regions occupied 
by the clouds in the remnant, and more than the flux den- 
sity of the entire source at that frequency. We conclude, 
then, that 7r°- decay is the dominant mechanism for 7-ray 
production in the new remnant. When an expanding SNR 
collides with a cloud, a fraction of the nuclear CR com- 
ponent is transported into the cloud by convection. This 
penetration produces a significant enhancement of the 7- 
ray luminosity because of the Tr'^-decays resulting from the 
p — p interactions. The 7-ray flux expected at the Earth 
is: 

F{E > 100 MeV) = -^-g^Mdoud d-^m^\ (2) 
47r 

where q-y is the 7-ray emissivity per H atom of the cloud 
of mass Mcioudj d is the distance to the source, and toh is 
the H-atom mass. 

Using the observed 7-ray fluxes from the CI and C2 
positions, the adopted distance of 470 pc, and our cloud 
mass estimates derived from the Hl-observations, we find 
that there is an enhancement in the CR energy density 
in the clouds with respect to the value observed near the 
Earth (Dermer 1986). The enhancement factor is calcu- 
lated to be ^ 4.5 for cloud CI and ^ 13 for cloud C2, 
clearly indicating that protons are being accelerated at 
the remnant shock front. 

Since at high energies the spectrum of Tr^-decay 7- 
rays mimics the spectrum of the parent protons, we can 
infer from EGRET observations that at the position of 
C2 the proton spectrum is N{Ep) oc ^-2-6±o.2 rj,^jg 
power-law index is similar to that of the electron com- 
ponent at the same position, as deduced from the radio 
brightness spectral index. The 7-ray spectrum towards CI 
may be harder: N{Ep) oc i<;-2-3±o.45^ r^j^g radio emission 
from this region is too weak to provide a reliable electron 
spectral index for comparison. In both cases the proton 
spectral index is consistent with the theoretical value of 
r = 2.42 that is expected from diffusive shock acceleration 
in shell- like remnants (Biermann & Strom 1993; see also 
Gaisser et al. 1998, who consider a harder injection spec- 
trum with r ~ 2.3). Leakage from the galactic disk results 
in the steeper i?^^-^^-spectrum observed near the Earth 
(Biermann & Strom 1993). We estimate that the proton 
high-energy cutoff, considering the age of the SNR, should 
be about 80 TeV (see Biermann & Strom 1993 for the for- 
mulae). 



5. Conclusions 

We have discovered a new southern SNR in the 
Capricornus region. It is a large and faint shell-type non- 
thermal radio source which has till now been obscured by 
confusion with the diffuse galactic radiation. There are 
three 7-ray sources in the 3EG catalog superimposed on 
this remnant. Two of these high-energy sources coincide 
with the positions of HI clouds detected at 21 cm in nar- 
row velocity intervals. We have estimated the masses of 
the clouds and calculated the CR flux enhancement at 
the SNR, showing convincing evidence of proton accel- 
eration at the shock front, as predicted by the theory. 
Interestingly, it is possible that the proton spectrum slope 
is not the same at all points in the remnant. The electron 
and proton components of the CR flux seem to display 
similar spectral slopes towards the single position where 
both spectra can be determined. Although the possibility 
of a pure chance association cannot be completely ruled 
out, it is clear that this southern cluster of 7-ray sources 
deserves more attention as a potential natural laboratory 
for CR studies. In particular, future X-ray and TeV ob- 
servations of the new Capricornus SNR are required to 
provide the information necessary to determine the up- 
per cutoffs in the spectrum of both the leptonic and the 
hadronic CR components within the source. 
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